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1,0  INTRODUCTION 

During  the  first  three  or  four  months  of  CY  1976,  the 
Defense  Nuclear  Agency  organized  and  funded  an  extensive  program 
of  one-dimensional  and  two-dimensional  computer  code  calculations 
dealing  with  low-yield,  shallow-buried  nuclear  device  detonations. 
The  main  thrust  of  this  effort  was  to  better  define  the  blast, 
thermal  and  initial  nuclear  radiation  environments  associated 
with  such  detonations.  A secondary  objective  was  to  attempt  to 
compare  the  results  obtained  from  several  different  complex  2-D 
rad-hydro  codes  which  are  being  run  around  the  country. 

In  order  to  make  these  comparisons  of  code  results  meaning- 
ful, DMA  requested  that  ERDA-Livermore  generate  the  source  tapes 
to  be  used  by  all  the  participants.  Responding  to  that  request, 
Livermore  (Jeff  Thomsen)  supplied  source  tapes  for  each  depth  of 
burial  considered. 

As  a portion  of  a Kaman  Sciences'  investigation  for  DNA  of 
the  test  design  for  an  underground  cavity  experiment  (HURON 
KING) , it  became  necessary  for  Kaman  to  perform  2-D  rad-hydro 
computations  for  two  of  the  geometries  of  interest.  The 
computations  were  done  using  the  KSC  version  of  the  RAD9  Code. 

This  report  presents  some  of  the  results  of  those  two  calcula- 
tions . 

Following  this  introduction,  the  RAD9  Code  will  be  described 
briefly.  Then,  the  principal  results  will  be  presented  and 
discussed;  for  individuals  who  are  interested  in  the  detailed 
results,  more  data  are  included  in  the  Appendix.  The  report 
concludes  with  some  suggestions  for  further  work,  using  the  RAD9 
Code . 
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2,0  DESCRIPTION  OF  THE  RAD9  CODE 


The  Kaman  Sciences'  RAD9  Code  is  a two-dimensional  multi- 
material hydrocode  which  treats  material  strength  with  elastic- 
plastic  relations  and  handles  radiation  flow  (and  conduction) 
effects  with  a single-group  greybody  formulation.  Hydrodynamic 
difference  equations  are  explicit  and  accurate  to  terms  of 
first  order,  while  radiation  is  generally  implicit.  Massless 
tracer  particles  are  used  to  define  surface  positions  and  move 
across  the  Eulerian  grid,  so  that  the  code  yields  a Lagrangian- 
type  definition  of  moving  surfaces.  The  code  is  capable  of 
handling  up  to  nine  materials  with  two  materials  per  cell. 
Transfer  of  material  between  mixed  Eulerian  cells  is  effected 
by  forcing  a common  velocity  (and  pressure)  on  the  two 
constituents,  computing  the  total  increase  in  internal  energy 
within  the  cell,  and  partitioning  this  energy  between  the 
materials  by  assuring  that  energy  increments  are  proportional 
to  fractional  volumes  occupied  by  the  various  materials. 

These  fractional  volumes  are  computed  by  dividing  the  mass  of 
each  material  by  its  density. 

Generally,  Tillotson's  equations  of  state  are  used  in  the 
code;  however,  in  some  calculations  other  analytic  or  tabular 
forms  are  used.  A Erode  fit  to  the  Hillsenrath  data  has  been 
used  for  air  and  tabular  data  from  Sandia  Corporation  (CHART-D) 
has  been  employed  for  various  materials,  particularly  those 
used  in  nuclear  weapons.  Greybody  opacities  used  in  the  code 
are  obtained  from  Heubner's  (LASL)  tabular  data  and  they  include 
line  spectra. 
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3.0 


RAD9  CODE  RESULTS 


Some  selected  results  from  the  KSC  RAD9  Code  runs  for  the 
0.5-meter  and  3-meter  depth  of  burial,  1-KT  source,  are 
presented  in  this  section.  First,  the  maps  of  the  motion  of 
the  various  materials  in  the  problem  are  shown  for  specific 
times  after  the  burst.  Next,  some  temperature  profiles  are 
shown  for  the  3-meter  calculation  only,  and,  finally,  the 
partition  of  energy  between  the  soil  and  the  air  is  presented 
in  tabular  form. 

3.1  Material  Profiles 

The  material  profiles  corresponding  to  the  0.5-meter  DOB 

lation  are  shown  in  Figures  3.1  through  3.3.  The  first 
i--rofile  at  t = 0.175  msec  (Figure  3.1)  shows  the  characteristic 
shape  of  the  early-time  material  "bubble".  It  should  be  noted 
that  although  the  bubble  is  truly  oblong  in  shape  (larger 
vertically  than  horizontally) , the  distortion  in  shape  is 
somewhat  accentuated  by  the  choice  of  scales  on  the  plots. 
Reference  to  Figures  3.2  and  3.3,  which  show  material  profiles 
representative  of  later  times,  indicates  that  the  profile  evolves 
as  if  the  source  were  expanding  from  a region  '!  to  6 meters  above 
the  initial  soil/air  interface.  It  is  also  apparent  that  a large 
volume  of  soil  is  moved  upward  and  outward  at  these  relatively 
early  times,  supporting  the  expanding  air  blast  pressure  wave. 
Note  that  near  the  soil/air  interface  there  is  a cusp,  which 
is  again  characteristic  of  a low  air  burst  geometry  rather  than 
a buried  burst  geometry. 

The  material  profiles  corresponding  to  the  3-meter  DOB 
calculation  are  shown  in  Figures  3.4  through  3.7.  As  one 
would  expect,  the  profiles  shown  here  for  the  deeper  buried 
burst  are  significantly  different  from  those  shown  for  the 
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FIGURE  3.1  RAD9  CALCULATION  MATERIAL  PROFILES,  1-KT  AT 
0.5-METER  DEPTH,  t = 0 AND  t = 0.175  MSEC 
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FIGURE  3.2  RAD9  CALCULATION  MATERIAL  PROFILES,  1-KT  AT 

0.5-METER  DEPTH,  t = 0.180  MSEC  AND  t = 0.193 
MSEC 
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0.5-meter  case.  The  early-time  expansion  (Figures  3.4  and  3.5) 
is  characteristic  of  the  so-called  "earth-piston"  effect  which 
has  been  observed  often  when  large  high-explosive  (H.E.) 
charges  are  buried  and  detonated.  For  the  later  times 
(Figures  3.6  and  3.7),  the  piston  appears  to  push  up  through 
the  top  and  expand  more  significantly  in  the  vertical  direction. 
To  be  specific,  at  t = 1,734  msec  (Figure  3.5),  the  piston 
dimension  along  the  ground  surface  is  about  8 meters  and  close 
to  8 meters  above  ground  zero;  however,  at  t = 3.965  msec 
(Figure  3.7),  the  piston  size  along  the  ground  has  expanded  to 
about  12  meters  whereas  its  height  is  more  than  19  meters. 

When  one  refers  to  the  air  blast  pressure  contours  (Appendix) 
it  is  apparent  that  this  expansion  of  the  soil  piston  has  a 
significant  effect  upon  the  pressures. 

3.2  Temperature  Profiles 

The  temperature  profiles  at  selected  times  after  burst  for 
the  3-meter  depth  calculation  only  are  included  in  Figures  3.8 
and  3.9.  The  temperatures  are  given  in  eV.  units  where  1 eV . 
is  about  11,800°K,  In  each  case,  the  profile  defines  the 
maximum  space  extent  of  the  labeled  temperature  given;  that  is, 
all  temperatures  inside  the  profile  are  higher  than  the  label 
value . 

The  general  behavior  illustrated  in  the  figures  shows  a 
very  hot  region  of  soil/iron  mixture  at  early  times  (t  = 1.297 
msec  in  Figure  3.8)  located  slightly  above  the  air/soil  inter- 
face. As  time  passes,  this  hot  region  cools  somewhat,  expands, 
and  moves  upward  - until  at  t = 3.965  msec  (Figure  3.9),  it 
extends  over  a region  of  about  12  meters  vertically  which  is 
centered  near  a height  of  12  meters  above  ground  zero.  The 
temperature  of  this  region  is  about  3 eV,  (or  about  30,000°K) 
near  its  outer  surface  and  hotter  on  the  interior.  The 
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FIGURE  3.4  RAD9  CALCULATION  MATERIAL  PROFILES,  1-KT  AT 
3-MErrER  DEPTH,  t = 0 and  t =1,297  MSEC 
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FIGURK  3.5  RAD9  CALCUI.ATTON  MATERTAI-  FROFILFF,  1-KT  AT 


3-MCTKR  DEPTH,  t = 1.406  MSEC  AND  t = 1.734 
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figure  3.8  RAD9  CALCULATION  TEMPERATURE  PROFILES,  1-KT  AT 
3-METER  DEPTH,  t = 1.297  MSEC  AND  t = 1.93 


MHEKS 

FIGURE  3.9  RAD9  CALCULATION  TEMPERATURE  PROFILE 


AT  3-METER  DEPTH,  t 


3.965  MSEC 


general  behavior  illustrated  in  Figures  3.8  and  3.9  indicates 
that,  concerning  thermal  radiation  heating  of  the  air  and  the 
ground  surface  at  later  times,  the  3-meter  depth  burst  acts 
more  like  an  air  burst  at  about  12-meter  HOB  and  reduced  yield. 
The  thermal  yield  would  be  reduced  due  to  the  large  amount  of 
soil  moved  up  and  out  which  provides  an  effective  shield  to 
the  radiation  emitted  by  the  hot  fireball  which  is  characterized 
by  the  temperature  profiles.  This  shielding  phenomenon  adds 
complexity  to  the  already-dif f icult  problem  of  computing  the 
total  thermal  radiation  (vs.  time)  reaching  the  ground  surface 
at  various  ranges  from  GZ . 

3.3  Energy  and  Mass  Results 

Tables  3.1  and  3.2  include  summaries  of  the  energy  and 
mass  results  obtained  from  the  0.5-meter  and  3-meter  depth  of 
burst  calculations,  respectively.  For  the  0.5-meter  case,  the 
ERDA-Livermore  source  tape  indicated  that,  at  the  time  the 
RAD9  calculation  was  started  (0.175  msec)  the  contribution  from 
the  iron  to  the  total  mass  and  energy  is  only  about  0.06%  and 
decreases  with  time;  for  this  reason  the  iron  data  were  omitted 
from  Table  3.1. 

Reference  to  Table  3.1  indicates  that  the  fraction  of  the 
total  energy  going  into  the  air  is  increasing  with  time. 

However,  it  is  also  evident  that  this  calculation  would  have 
to  be  carried  out  to  times  of  about  10  msec  to  accurately 
determine  the  air/soil  energy  partition.  The  table  listings 
indicate  that  energy  and  mass  is  conserved  very  well  during 
the  calculation. 

Table  3.2  listings  show  that  the  fraction  of  total  energy 
partitioned  to  the  air  for  the  3-meter  buried  charge  is  much 
less  than  for  the  0.5-meter  case,  which  is  as  expected.  Also, 
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TABLE  3.1 


SUMMARY  OF  ENERGY  AND  MASS  RESULTS  FROM  0.5-METER 
DOB  KAMAN  SCIENCES  RAD9  CALCULATION 


MATERIAL 


Soil 

Air 

Totals 


Soil 

Air 

Totals 


Soil 

Air 

Totals 


Soil 

Air 

Totals 


Soil 

Air 

Totals 


Soil 

Air 

Totals 


Soil 

Air 

Totals 


INT.  E 
(jerks) 


KIN.  E 
(jerks) 


^jerk  = 10^^  ergs 


*1.579(+10)  = 1.579x10^® 


TOTAL  E 
(jerks) 


TOTAL  MASS 
(gms) 

TIME 

(msec) 

1.579  ( + 10)*: 
1.99 (+7) 
1.581  ( + 10) 

0.175 

1.579 (+10) 
1.99 (+7) 
1.581 (+10) 

0.180 

1.579 (+10) 
1.99  (+7) 
1.581 (+10) 

0.193 

1.579 (+10) 
1.99 (+7) 
1.581  (+10) 

0.248 

1.579  ( + 10) 
1. 99 (+7) 
1.581  (+10) 

0.342 

1.579 (+10) 
1.98  (+7) 
1.581  ( + 10) 

0.469 

1.579  ( + 10) 
1.98  (+7) 
1.581  (+10) 

0.472 

TABLE  3.2 


SUMMARY  OF  ENERGY  AND  MASS  RESULTS  FROM  3-METER  DOB 
KAMAN  SCIENCES  RAD9  CALCULATION 


MATERIAL 


INT.  E 
(jerks) 


KIN.  E 
(jerks) 


0.13 

1869 

13.22 

1883 


0.078 

1811 

16.30 

1827 


TOTAL  E 
(jerks ) 


15.0 

4575 

18.8 

4608 


5.39 

4585 

17.68 

4608 


0.46 

4538 

62.68 

4601 


.31 

496 

.3 

574 


0.31 

4495 

77.72 

4573 


TOTAL  MASS 
(gms) 


1.76 (+4) 
2.535  (+10) 
2.32  (+7) 
2.537  ( + 10) 


2.55  (+4) 
2.535 (+10) 
2.32  (+7) 
2.537  ( + 10) 


2.01 (+4) 
2.535  (+10) 
2.32  (+7) 
2.537  (+10) 


1.53  (+4) 
2.5346 (+10) 
2.32  (+7) 
2.537  (+10) 


TIME 

(msec) 


1.297 


^jerk  = 10^®  ergs 
*1.76  (+4)  = 1.76x10'* 
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it  is  evident  that,  at  early  times,  the  contribution  of  the 
energy  from  the  iron  to  the  total  energy  is  significant;  as  time 
proceeds,  this  contribution  rapidly  decreases.  It  is  estimated 
that  this  calculation  would  need  to  be  carried  out  to  times  near 
100  msec  to  obtain  accurate  data  regarding  air/soil  energy 
partition.  Using  the  data  available  from  these  calculations, 
an  "educated  guess"  would  say  that  the  fraction  of  energy 
appearing  as  air  blast  for  the  3-meter  case  is  less  than  one 
would  predict  using  the  standard  effective  yield  concept.  This 
conclusion  should  be  verified  by  carrying  out  the  calculations 
to  the  required  later  times. 

Carrying  out  the  KSC  RAD9  calculations  of  the  0.5-meter  and 
3-meter  cases  to  later  times  would  not  require  excessive  computer 
costs.  The  results  reported  here  were  obtained  using  a CDC  7600 
Computer  with  the  following  run  times: 

(a)  0.5-meter  case  - 19  min. 

(b)  3-meter  case  - 16  min. 

3.4  Contour  Plots  of  the  Results 

Plots  of  space  contours  of  pressure,  density,  velocity 
vectors  and  internal  energy  are  presented  in  an  Appendix  to 
this  report.  The  plots  are  given  for: 

1,  0.5-meter  DOB  (1-KT) ; t = 0.472  msec 

2.  3-meter  DOB  (1-KT);  t = 3.981  msec. 

The  Z-  and  R-axes  are  labeled  in  centimeters;  for  the 
0.5-meter  case  the  air/soil  interface  is  at  Z = 1000  centimeters 
and  for  the  3-meter  case  the  air/soil  interface  is  at  Z = 1600 
centimeters . 

For  the  0.5-meter  case,  some  pressure  vs.  radius  plots 
are  presented,  corresponding  to  selected  values  of  constant  Z. 
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4.0  SUGGESTED  FURTHER  WORK 

It  is  suggested  that  the  following  work  be  done  with  the 
KSC  RAD9  Code: 

1.  Continue  the  KSC  RAD9  calculations  of  the  1 KT- 
3 meter  case  to  about  200  msec  and  the  1 KT-0.5 
meter  case  to  about  20  msec  to  determine  the 
yield  fraction  contributing  to  air  blast  in  each 
case . 

2.  Using  the  results  obtained  from  1.  above,  compare 
these  results  with  the  AFWL  HULL  Code  calculations 
(using  the  Tillotson  Equation  of  State)  for  the 
same  cases. 

3.  Using  the  KSC  RAD9  Code,  conplete  the  following 
sensitivity  studies; 

(a)  Determine  the  effect  of  changing  (within 
sensible  limits)  the  linear  and/or  quadratic 
terms  in  the  artificial  viscosity  used  in 
the  problem.  Most  of  the  codes  used  for 
these  calculations  incorporate  the  artificial 
viscosity  formalism. 

(b)  Determine  the  effect  of  zone  size  near  the 
source  and  at  the  shock  front.  There  is  some 
evidence  that  even  the  "shock-following" 
routine  used  by  AFWL  is  not  adequate  to  define 
the  peak  overpressure  accurately  at  low  (less 
than  20  psi)  overpressures.  Kaman  Sciences  is 
developing  a "packaging"  routine  which  will 
help  to  define  these  peak  overpressures  without 
the  use  of  excessive  computer  time. 
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(c)  Determine  the  parameters  which  cause  the  most 
severe  radiation  leakage  in  the  calculations 
considered.  One  candidate  is  the  value 
chosen  for  the  "flux-limiter"  variable. 

(d)  Determine  the  effect  of  allowing  materials  to 
"mix"  inside  the  Euler  cells,  particularly 
when  high  opacity  materials  are  adjacent  to 
low  opacity  materials.  This  situation  could 

lead  to  excess  leakage  of  radiation,  significantly 
influencing  the  results. 
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APPENDIX 


PLOTS  OF  THE  RAD9  CODE  RESULTS 


This  appendix  includes  Figures  A.l  through  A. 40.  The 
subjects  covered  can  be  summarized  as: 

1 . 1-KT;  0.5-meter  DOB;  t = 0 .472  msec 


(a) 

Figures  A.l 

through  A. 7: 

Pressure  Contours 

(Z-R  space) 

(b) 

Figures  A. 8 

through  A. 10 : 

Pressure 

vs . R 

Profiles 

(c) 

Figures  A. 11 

through  A. 17: 

Density 

Contours 

(Z-R  space) 

(d) 

Figure  A. 18: 

Velocity  Vectors  (Z-R 

space) 

(e) 

Figures  A. 19 

through  A. 22: 

Internal 

Energy 

Contours  (Z- 

R space) . 

1-KT 

; 3-meter  DOB 

; t = 3.981  msec 

(a) 

Figures  A. 23 

through  A. 29: 

Pressure 

Contours 

(Z-R  space) 

(b) 

Figures  A. 30 

through  A. 36: 

Density 

Contours 

(Z-R  space) 

(c) 

Figure  A. 37 : 

Velocity  Vectors  (Z-R 

space) 

(d) 

Figures  A . 38 

(Z-R  space) . 

through  A. 40: 

Internal 

Energy 
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.472  MSEC;  CONTOURS  OF  CONSTANT  PRESSURE 


FIGURE  A. 2 RAD9  CALCULATION,  1-KT  AT  0.5-METER  DEPTH, 
t = 0.472  MSEC?  14.7  PSI  PRESSURE  CONTOUR 
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figure  a. 12  RAD9  CALCULATION,  1-KT  AT  0.5-METER  DEPTH, 
t = 0.472  MSEC;  1x10  g/cc  DENSITY  CONTOUR 
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FIGURE  A. 16  RAD9  CALCULATION,  1-KT  AT  0.5-METER  DEPTH 
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CONTOUR  lEn  LINE'i  SLIDE 


t = 0.472  MSEC;  1x10  ERGS/GM  INTERNAL 
ENERGY  CONTOUR 
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.981  MSEC;  CONTOURS  OF  CONSTANT 
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FIGURE  A. 34  RAD9  CALCULATION,  1-KT  AT  3-METER  DEPTH, 
t = 3.981  MSEC;  3xl0~^  g/cc  DENSITY 
CONTOUR 
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FIGURE  A. 36  RAD  9 CALCULATION,  1-KT  AT  3-METER  DEPTH 
t = 3,981  MSEC;  1.0  g/cc  DENSITY  CONTOUR 
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figure  a. 38  RAD9  CALCULATION,  1— KT  AT  3— METER  DEPTH, 
t = 3.981  MSEC;  CONTOURS  OF  CONSTANT 
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Nevada  Operat  Ions  Office 

ATTN:  Doc.  Ci*n.  tor  Tecii.  I-ib. 

I'nion  Carbide  Corporation 
Holitleld  Natlon.il  L.iboratory 

.ATTN:  Do.  . Con.  for  Tech.  Lib. 

ATTN:  Civil  Def.  Res.  Pro.). 

OTHER  (iOVERNMKNl  A(;KNCIKS 

Dep.irtment  ol  the  Interior 
Bureau  of  Mines 

ATTN:  Tech.  Lib. 

Department  ol  the  Interior 
I’S  (a'ological  Survey 

ATTN:  Cecil  B.  Raleigh 
ATTN:  .1.  H.  Healy 

DEPARTMENT  OF  DEFENSE  pONTRACTORS 

Aerospiice  Corpi»ral  Ion 

ATTN:  Tech.  Info.  Services 

Agbablan  Associates 

ATTN:  M.  Agbabian 

Applied  Theory,  Inc. 

2 cy  ATTN:  .lohn  G.  Trnllo 

Avco  Research  ft  Systems  (Jroiip 

ATTN:  Rese.\rch  Lib..  AH  JO,  Km.  7201 

Battelle  Memorial  Institute 

AITN:  Te.hntcal  l.lhrary 

The  BDM  la>rpi)rat  Ion 

ATTN:  Techrjical  Llbrarv 


DEPARTMENT  J\F  DEFENSE  C0N'i%\CT0RS^  i 

The  B«>eing  IT'mpany 

ATTN:  Aerospace  Library 

ATTN:  R.  M.  Schmidt 

l!.aIiiornia  Research  & Technology,  Inc. 

ATTN:  Sheldon  Shuster 

ATTN:  Ken  Kreyenhagen 
ATTN:  Technical  Library 

Calspan  Corpor.u  ion 

ATTN:  Technical  Library 

C I V 1 1 /Nuc lear  Systems  Corp. 

ATTN:  Robert  ('rawt(>rd 

Pniversiiy  of  Dayton 

Indijstrial  Security  Super  KL-b05 
ATTN;  Hal  lock  F.  Swilt 

I'nlversitv  ol  Denver 

(ailorado  Seminary 

Denver  Research  Institute 

ATTN:  Sec.  Oflicer  for  d.  Wisi>tski 

F.<>t»(>,  Inc. 

Albuquerque  Division 

ATTN:-  Technical  Library 

Card,  I nct»rporat  ed 

ATTN:  G.  L.  Neldh.irdt 

Genera  I Klee  I r Ic  Compiuiy 

TEMi’O-Cent  er  for  Advanced  Snjdies 
ATTN;  DAS  I AC 

I IT  Rese.trch  Inst  lt\Jle 

ATTN:  Technical  Library 

Institute  tor  Detense  Analyses 

ATTN:  IDA  Librarian,  Ruth  S.  .Smith 

Kaman  AviDvne 

Dlvi.sion  ol  Kaman  Sciences  Corp. 

ATTN:  F.  S.  ('rlscione 
ATTN:  Technical  Library 

Kaman  Sciences  Corp.'r.il  ion 
ATTN:  Library 
ATTN:  Don.ild  C.  Sacljs 

Lockheed  MIssili*s  A Space  Companv,  Inc. 

ATTN:  Te.hnlcal  library 

I.ockiieed  Missili-s  ^ Space  (a>mpany»  Inc. 

ATTN:  Tom  Geers,  D/')2-13,  Bldg.  205 
ATTN:  Tech.  Info.  Center  D/CGLI. 

McDonnell  Dcniglas  Corporal  ion 
ATTN:  Robert  W.  Halprln 

M.'rrllt  CASES,  Incorp.'ial.'d 

ATTN;  Te.hnlcal  Llbrarv 
ATTN:  .1.  L.  Merritt 

The  Mitre  Corporation 
ATTN:  library 


67 


1 


DKlVXR’fMKNT 

OF  DKFKNSK  CONTR^UTORS  (Continuod) 

DF.PARTMKNT  OF  DKFKNSK  CONTRACTORS  (Continued) 

N.uhan  M.  NewTuirk 

TRW  Jietonso  i Spat  e Sy.s.  Crtiup 

Consul  t in^; 

Knginoorln^;  Sorvloos 

ATTN:  T.'ch.  lull).  (Vnl I'l/S- 1 9 10 

BI06A  Civil 

Kn>»inoorinK  Building 

ArrN:  D.  H.  B.ior.  Rl-2nh 

Cniversitv 

ol  111  intMs 

ATTN:  R.  K.  Plobuoh,  R1-J07H 

ATTN : 

Nathan  M.  Nowni.jrk 

ATTN:  1.  K.  Alber.  Rl-IOOH 

? ly  ATTN:  I’.tiT  K.  Il.ii,  Kl/2170 

Physios  Intoriuit  iotuil  Ci>mp.iny 

ATTN: 

t)t»o . Ci'u.  for  Tooh.  I.ib. 

rkW  Defense  it  Sp.Jt'o  Sys.  Croup 

A r I'N : 

Doc.  Ci>n.  It'r  Rohort  Swilt 

San  Bern.trdlno  i)por.it  ions 

ATTN: 

Doo.  Ci»n.  tor  Donnis  Orph.il 

ATTN:  K.  V.  WoilR,  927/7K’ 

ATTN 

Dt'o.  Con.  for  L.irrv  A.  Bohrmatni 

A ITN : 

Dot'.  Ct*n.  tor  Frod  M.  SaiuM- 

Cn Ivor sal  Analytios.  Ino. 

ATTN  ■ 

Dtu*.  Ct»n.  tor  Charlos  Ctidlroy 

ATTN:  i:.  I.  Field 

ATTN : 

l)t»o . Ci’n.  ft>r  F.  T.  M«’oro 

CKS  Rosearrh  Company 

K & n Assoi 

i.ilos 

ATTN:  loohnioal  Libr.iry 

ATTN : 

U'itlian)  B.  Wrigiit,  Jr. 

ATTN : 

J.  {•.  Low  is 

The  Kri  H.  W.ing 

ATTN : 

Toohnioa 1 Library 

I’ivil  K.'  inooring  R.soli.  Kao. 

ATTN 

Rt>hort  Tort 

I'niversi  t ■ St.it  ion 

A'l'lN 

Hon>  V Ci>t»por 

ATTN:  No.il  Baum 

ATTN 

H.irtTil  1..  Brodo 

ATTN:  Larry  Hloklo 

ATTN 

Cvrus  P.  Knowles 

A ITN 

.lorrv  Carponlor 

W.ishingit'n  St. tie  CnIvorsity 

A TTN 

Jorrv  Stookti»n 

Ailministrat  ivo  ('rganizat  Itui 

ATTN:  Arthur  Milos  H«>liort  for  Ceorgo  Duval 

So  liMU’i'  AppI  io.u  ions,  Ino. 

ATTN 

n.  F.  M.ixw.'M 

Woidlingor  Assoo.  (’tinsuiting  Enginoors 

ATTN 

Daviil  Bornsloln 

ATTN:  .1.  W.  Wright 

A'lTN:  Melvin  L.  Raton 

Soionoo  Appl  io.u  i.>ns.  Ino. 

ATTN 

Tot  hnli  al  Mbrarv 

Weidlingor  Asst»o.  Ct'nsuItIng  Kngineers 

ATTN:  .1.  l.-,onberg 

St'uthwost  Rt'so.in  h Insiituto 

ATTN 

Will  rod  K.  Baker 

Wt‘st  inghouse  Fleet  t io  Corp. 

AITN 

A.  H.  Win/el 

M.ir  ino  Division 

AITN:  W.  A.  Voir 

St.inlord  Rosoaroh  Insiituto 

ATTN 

Burt  R.  Cast  on 

ATTN 

Ceorgo  K.  Abrahamst>n 

Sv.stoms,  St 

lonoo  StOtw.iro,  Ino. 

ATTN 

Tluunas  1).  Kinov 

ATTN 

Tod  ('horrv 

AT  I'N 

Toohnioal  Library 

Al’TN 

Dtin.iltl  R.  Crine 

Toria  Tok. 

Ino . 

A TTN 

Sitlnoy  Croon 

ATTN 

Toohnioal  Library 

Toti.i  Toth 

, Iiu  . 

ATTN 

Toohnioal  Library 

ATTN 

LI  -San  ilw.ing 

68 


